gH Dike eres = went bie 


mit babi Dan 


SO Aiecre ine 


ot te tien nasheed 
Sheree 


2 lemon 
x Se Sa 


Shatin pak 


eirees 


ST ree 


hlbchalecliedaenapirnad tote te ete 


In Moone ran thn je oe Sw fe 


mjtaiinaie Nae Cerin 
beatin Titi kan 5 


—tsteliae Semi rone 


eee Uelcch oty 


Heerlen tne: Raina 


: rachael ee ee 


2S ee ae Nene 


im be a iw erie Se 


ha ee lh em hi 


2 Fw ne te Naeem tial nyt tape n pi 
hal he iia es 


taitoereintib-iee itsbntoierse ty 
Pe Fea 


baiehe he ima any. 


vl ~ stim Weetn Seilem esiigu ke 


esl = nee Nan ng 
ees pmete ate 
— Hz iten dae etd 


nite Sel “i aie eid os 
Stinbwcie Minti 


needed onda 


Cate etlge ein Neng 


lode Lar 
andere 


Sein She ionialail 


9 re mtewibew lon tie g, 


“aethoaliond slim Sa 


edie ol ane ae 


Persie) Me eis wee! 


eholeHwietia ca nwdenhose 


Cait i 


Prien eon 


fatale dw a6 ature eioe 
Weleederucatietanercen ee ere 
Rite petepere aden cy ke 


eel se rine VS ele 


erin ee 


folie sie entak emer 


Far 


Petites cate 


te Felt terebe wairana 


re H tonal ony 
TOA SPS eke le Ng 


neste nl wg te 


eee 


pereetl eee 


PL RR Sh relicricertens 


mS ye latinlhigs. Sanh 


Gnlion the ioe in gofienes 


Petes 


i eis 
ENE wae beetb of 


Sh tir eee Bas 
Rm niBenelte here, 


onan beens da pty 
Aim imine ie 2 


; aoeche 
inl nWenele Sw aac 
ie 


ellen eatretaneesiee” 
Hh esieen Hd Go's 


ln a 
" Sl rm Memety 
Sahih ic eae ete he teeter 


pale 


SP nbn od Pestle 


Sai 


barvicar crater 


re niiedieih othe wene 
MHL N Obie 2 Fi 


PIN a ers rel ielerart cain t pi 
eee 


What lime 


$e Metron g 


Gx wpnis 
UNIERSTTATIS 
AABERTAEASTS 


ag loa 


io nee 
tS Ae bi 


' 
Ves f 


THE UNIVERSITY OF ALBERTA 


CONVERGING FLOW 


OF A 
VISCOELASTIC FLUID 


by 


HIN-SUM LAW 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 
OF MASTER OF SCIENCE 
IN 
CHEMICAL ENGINEERING 


DEPARTMENT OF CHEMICAL ENGINEERING 
EDMONTON, ALBERTA 


SPRING, 1977 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


httos://archive.org/details/Law197/ 


ABSTRACT 


The converging flow patterns obtained when a 0.2% aqueous 
solution of Separan AP-30 flows from a large reservoir into capillary 
tubes of 0.2 and 0.3 cm diameter were studied by using a flow vizual- 
ization technique. Detailed observations of the variation of converging 
Flow patterns with fiow rate and tube diameter were made. In addition 
complete velocity profiles were determined quantitatively at several 
positicns upstream of the entrance. 

In common with a number of earlier studies this work shows that 
the stable converging flow fieid is in the shape of a wine-glass stem 
with circulating flows developed in the corner regions. Comparison of 
the velocity profiles and flow patterns with available models show these 
to be largely inadequate. The initial cone semi-angle (90 was found to 
be related to the inside diameter of the capillary tube and average fiow 


. velocity through the equation: 


ne 0.60 
Ps = 64.05 (7) 


The centerline axial velocities within the central core could be 


expressed as: 


His 0.545 aio 
aes = 1.645 exp (-4.04 ‘T) ip! 


Contrary to earlier studies, the velocity profiles showed 
considerable development prior to the contraction. At the contracticn 
velocity at the centerline was 95-100% of the centerline velocity for 


fully developed tube fiow. 
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CHAPTER I y 
INTRODUCTION 


bel General 
The converging flow of viscoelastic fluid is of interest in 
petroleum recovery schemes. In secondary recovery operations, fluid 
injection is commonly used as an external energy supplying the reservoir 
in ordér to recover the oil trapped in the porous rock. This external 
energy is needed because the initial natural energy in the reservoir such 
as the gravitational energy due to the difference in hydrostatic level 
Of producing well and-upper part of the reservoir will sooner or later be 
depleted. Polymer flooding is a fairly new fluid displacement technique, 
which has been shown to be more effective than water flooding to develop 
a stable high oil saturation zone simply because the polymer solution has 
higher viscosity than water (52). The mechanism of polymer flooding is 
not fully understood but is probably related to the flow behaviour of 
polymer solution in porous media. Actually, the type of flow through 
porous media is very similar to the converging flow into a capillary tube. 
On the other hand, the converging flow of viscoelastic fluid 
is of interest in many polymer processing operations sucn as extrusion 
and fiber spinning. The converging flow is one of the extensional 
flows since the fluid elements are subjected to large stretching deform- 
ations while accelerating into the capillary tube. Extensionai flows 


can be divided into two types as follows: 
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(i) forced extension 
In flow of this type, the shape of the material being 
extended is forced by the operation being performed. For 
example, the calendering where the material is forced to 


conform to the shape of the rollers as it passes through them. 


(47) free extension 

In flow of this type the shape of the extending material 
is not imposed by the operations, and the material chooses 
its shape subject to internal constraints. |More commonly 
continucis operations involve this type of extension which 
the material being extended has a lateral interface with a 
low-viscosity substance or is in lateral contact with portions 
of the same material] which are not undergoing extension. For 
example, the first case is illustrated in Figure (I-la) for 
spinning where the low-viscosity substance is air, and the 
latter case is illustrated in Figure (I-1b) for extrusion of 


polymer solutions. 


Since the converging flow is in the type of free extension, the 
experimental results of the converging flow may be similar or have appli- 
cation to fiber spinning and extrusion. 

Almost all polymer extension flow problems of interest to the 
rheologist involve free extension rather than forced extension. The 
extensional rheology of polymeric fluids has been considered by many 
investigators such as Balakrishnan and Gordon (4), Balmer (7), Cogswell 
(13, 14), Denson (17), Kizior and Seyer (29), Metzner and Metzner (31), 
Murch (33) and Oliver (39) for various types of free extension. These 


studies include information on the extension viscosity which is defined 
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4 
as the axial or normal stress in fiow direction diyided by the extensional 
strain rate. The extension viscosity was reported to be as much as a few 
orders of magnitude greater than the shear viscosity of the polymeric 
Fluvdssi4, 315-39). 

In the past little attention has been devoted to quantitative 
measurements of the velocity profiles of the extensional flows and relating 
these velocities to the parameters such as flow rate, diameter of orifice 
or een eliary tube and axial location along the direction of flow. The 
major majority of investigations of extensional flow have been concerned 
with the onset and prediction of flow instabilities because in many polymer 
processing operations, the rate of production is limited by the onset of 
flow instabilities. The studies of melt fracture, which is the name 
given to a gross distortion or waviness of extrudate ceformations that are 
Shown by polymer soiutions and melts flow through capillaries or slits at 
output above a certain criticai value, were reported by many investigators 
such as Bagley (2), Ballenger et al (5), Barnett (8), Blyler and Hart (9), 
den Otter (15), Giesekus (18), Han (22), Hurlimann and Knappe (23), 
McIntire (30), Pearson and Pickup (40), Rama Murthy (43), Rothenberger, 
McCoy and Denn (44), Southern and Paul (45), Tomita and Shimbo (46), 
Tordella (47, 48), Vinogradov and Manin (50) and White (56). Flow 
visualization techniques were commonly used by.almost all of the investi- 
gators. In addition, Kase (28) theoretically studied the stability of 
melt spinning and Petrie and Denn (41) have published review of the 
instabilities in polymer processing. 

A technique of measuring the velocity profiles of molten poly- 
mers during laminar flow was first represented by den Otter (16). 


Velocity profiles were measured for the fiow of some polymeric melts 
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(low density polyethylene, two high density polyethylenes and a polydi- 
methyl siloxane) through a rectangular slit. The velocities were deter- 
mined by measuring the time needed for a particle to cover a distance of 
one millimeter with a stopwatch. The velocity profiles measured were 

in good agreement with those calculated from the flow curves made by 
assuming a simple power law relationship over the whole shear rate range. 
The velocities were found to be zero at the wall in all cases even during 
melt fracture for a low density polyethylene and a polydimethyl siloxane. 

Metzner, et al (32) studied the kinematics of a converging 
velocity field in flow from a large reservoir into a smail tube. The 
velocity profiles were measured photographically by using smal] air 
bubbles as tracers. More details of their study will be discussed in 
the next section. 

Kanel (26) had studied the behavior of a viscoelastic fluid 
in "fiber spinning" and "Fano flows". "Fiber spinning" is a technique 
where the material is forced through a die at low velocity and is wound 
up at high velocity at some position downstream from the die. "Fano 
flow" is a thin column of fluid created by immersing one end of a thin 
glass tube in a reservoir containing polymer solution while the other 
end is connected to an evacuated vessel. When the level of the fluid 
dropped well below the tube end, a thin column of fluid would continue to 
flow upward and through the tube. By taking photographs of the 
extensional fluid column and making diameter measurements, Kanel suggested 
that the extensional flows were characterized by a history parameter P 
(when P=1 it was defined as K-] model). When two flows have the same 
value of the parameter P, they were seen to have the same deformation 


history over time intervals during which they cover the same range of 
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6 
extension rate. Kenel's (26) data and Weinberger's (54) data which Kenel 
used, indicated primarily a quadratic dependence of velocity on distance 
into the fluid column in Fano and fiber spinning flows. However a linear 
relationship which indicated a homogeneous steady extension was also noted 
in several instances. The solutions they used were different concen- 
trations of Separan AP-30 in a 50% glycerine - 50% water mixture. In view 
of the generally high fluid viscosities and the absence of significant 
shearing due to the low viscosity of the surrounding air, Kanel has assumed 
a flat velocity profile in his study. 

Kase (27) analyzed the velocity field within a molten spinning 
thread by solving the equations of continuity and momentum for Newtonian 
fluids. In solving the equations, the visocity was assumed known. He 
obtained a solution for the axial velocity which was an exponential 
function of the axial location. With the same reason as mentioned above, 
Kase has also assumed a flat velocity profile in his study. 

Although the above investigators have considered the velocity 
profiles within the extensional flow, only one of them (32) has measured 
the velocity profiles directly. Thus the goal of this study is to 
measure the velocity field in converging flow from a large reservior into 
a small capillary tube. These results may have application to the 


petroleum recovery schemes and polymer processing operations. 


ied Description of Convercing Flows of Viscoelastic Fluids 

A considerable number of experimental investigations have 
reported on the qualitative behavior of viscoelastic. fluids such as molten 
polymers in capillary rheometers. Results of these investigations have 
indicated that at low flow rates, the materials show a full 180° conver- 


gence with no evidence of secondary flow as a Newtonian fluid does (6,19,43). 
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7 
In other words, the flow at any point in the reservior is directed toward 
the tube entrance. At higher flow rates, circulating flows develop in 
the corner regions, their size increases as flow rate increases (6, 24, 25, 
32, 43). This velocity field of the converging flow of viscoelastic 
fluid is shown in Figure (I-2). At very high flow rate, a breakdown of 
this entrance velocity field occurs. 

The sink flow model has been used previously to describe the 
velocity field within the central core for converging flow by Uebler (49) 
and Metzner, et al (32). In their work, the velocity profiles within 
the central flow core for 0.5% Separan AP-30 solution, which flowed from 
an 18" wide square reservoir into a 1.48 in I.D. tube, were measured with 
flow rates of 457, 915 and 1372 cm/sec. The extent of the central core 
as measured by $ in Figure (1-2) was 10° to 15° from the extended center- 
line of the tube. Within this region the streamlines were straight and 
if extrapolated downstream from the tube entrance, converged to a point 
on the centerline of tne tube. Thus the flow kinematics of a viscoelastic 
fluid in this geometry was represented as that of a radial flow towards a 
virtual point sink located downstream from the tube entrance on the center- 
line of the tube. Uebler's results showed that the radial component of 
the fluid velocity within the central core was dependent only upon radial 
position; in other words, no shearing was present in that region. 

The initial cone semi-angle ($0) in Figure (I-2}, which is 
defined as the angle formed by the line tangent to the outermost stream- 
line at the contraction, has been studied by Ballenger and White (6) and 
Oliver and Bragg (35, 37, 38). Ballenger and White have measured the 
initial cone semi-angle for the flow of low-density polyethylene, poly- 


styrene, and isotactic polypropylene at 180°C in the reservoir approach 
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Figure I-2 


VELOCITY FIELD OF THE CONVERGING FLOW OF VISCOELASTIC FLUID 


Les pillary Tube 


NNANNANUNNRTRAN 


SS 


ww LREES 
SS 


ENSSSS 


Fee ee py ae ee ee ee ees ere eee a 
S ee _ 
ey J 
m 


CHM 


4 


Viigtae  F 


5 ~ ie | * = 4 
~ Rete aa ett sey 8 x 
ey pe oe ep Diane omagrirny-eey oan aaa auch 
a Vay , i A Ae ’ 
; ‘ Pe , a ae. i A 
gg By de ff 7 Fy 
4 x 7 


, ioe ? ig le ee . - 
Ph oye if ad # Pag Pi x » F- at 
ie ae eee eet eee eS Oe 
a ‘ rs ; 
ye ' 
i. "hy eae Ee! Cis 
rast . ns * ie 
ee He Hy 
ry * ’ %\ - 
=. Page: as oe ‘es 
ke La ct 
‘e a. % 
+. '_o 
4 


ra 


, A 5 
a eaten 
% at 


: 4 
| Z. san rig oe ee 
- ~ \% *\ <j : 
te 
‘ % Ng = 
Big sh —_ ” 
L Bs BONS m St 
‘Mc we 
i 7 ‘%\ h *. ) 
en 
rT t. *% 7 * © 
i 
; \. r 
ee 
r ', ey hy. . 
- eNO \E eo 
; Ss Fh ; 
a. \ ~ ™ 
~~ < *S 


RE ¢ 
ONS a 
Mp. * | q 
, A 
/ aa Was | a 

* * ee 
La Nace ie 

“A 2 ae 

. hei 
ry . Ss vt 


ie 
i mAs} 
sit ime the SAT S3 
N : 1 
{ 
7 UJ 
Hi 
< aay 17 
bi , Al a ; 
A ht iy 


g 
to a capillary extrusion rheometer. Initial cone semi-angle was related 
to the Weissenberg number which is defined as the ratio of the entrance 
pressure drop to the capillary wall shear stress. Oliver (36, 37, 38) 
has measured the initial cone cemi-angle for three concentrations of 0.1%, 
0.03% and 0.10% ET 597 solutions. Initial cone semi-angle was related 
to the average velocity of flow through orifice and the radius of the 
orifice. Cone angles were studied in a 3" wide square reservoir fitted 
with different orifices of diameters 0.1, 0.2, 0.5, 0.95 and 1.27 cm, and 
with average velocities of flow through the orifices from 7 to 800 cm/sec. 

Even today the important aspect of this eeveraitiq velocity 
field such as the defpendence of axiai velocity on axial location is not 
fully understood. For example, the flat velocity profile assumed within 
the central core by Metzner, et al (32) is suspect because of the high 
viscosity of the circulating region which surrounds the central core. 
Therefore, a detailed quantitative study of the velocity field of the 
converging flow of viscoelastic fluid is of interest in this study. In 
the present study, the velocity profiles, the cone semi-angles, and some 
rheological properties such as wall shear stress and stretch rate were 
measured for the converging flow field. This flow field was obtained 
by a 0.2% Separan AP-30 solution flowed from : 2.64 cm wide square 
reservoir into different capillary tubes of I.D. 0.20 and 0.30 cm with 


three different flow rates of 0.1598, 0.3196 and 0.6392 cite y/Se0% 
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CHAPTER II 


THEORY 


II.1 Introduction 

As mentioned in Chapter I, only a few investigators such as 
Metzner, et al (32), Kanel (26) and Kase (27) had studied the velocity 
profiles within extensional flow. The sink flow model suggested by 
Metzner, et al is considered first to analyze the results in this study 
because the type of converging flow and polymer solutions used in their 
study were nearly the same as this study. The empirical K-1 model 
represented by Kanel’ for "Fano flow" and "fiber spinning", and the 
exponential decay model represented by Kase for melt spinning are also 
considered because the shape of "Fano flow" and "melt spinning" are very 
Similar to the shape of the central core in this study. Both of these 
are free extensions and have the diameter of the fluid column or core 
decreasing along the direction of flow. These three models will be 


discussed separately in the following sections. 


II.2 Sink Flow Model 

The velocity field within the central core can be analyzed employ- 
ing a spherical polar coordinate system with its origin located down- 
stream on the centerline of the capillary tube at a point determined by 
extending the converging streamlines to intersection. Distances upstream 
of this origin are denoted by r* and the diameter of the capillary tube 
by D. The coordinate systems used for sink flow model, in analysis of 


the velocity field are shown in Figure (II-1). 


Yo * 4 re 
yen 
Uae + p 
a ft Hl : mil 
one en ' 
er Ay 


a dive 2rd: 


x 


a fa ey oid vetbuse bed e 


Pa ~ 


rebashe eos ee + aa 

ti Pay = = va Ce i. bs 

7 ea Ot Teli Set ee. 

ae hae rae xe hee * bolt 5: a 3 at rin ; “ 

wince thee eo ae te 

cae arr egaes ee ree cre e fi o ‘ oe 
“aoe biG. INS Tarr is ve | 10 

“5 ore Bie - : 

ine = = ¥ 


z Be nh: b eae 
18 O96 eehonige 3) oa" Bem 
eae ' Nae f { ah 
cue a a 4 P a } 
#236 Ry Gate (eid ne Lae 
mas i) a os eee 4 me: 
dans ago temlren: | elt te ager ok evEd: é ban 
aur Ses: E tae ie) ¥ ew v ne : 
ve ig *) emo Fea, oe g ‘ 
ed fF bi abo Det ae = al 0 no Yop bin 
S : abe = ork ‘ 7 Ae . : 
eed Dry: , fae veal 07 | 
a eg 4 bY ie ae rea Seer via 
# a4 e Pa ff ~ 
¥ x - s io 
4 = +) 
a hea rt “al 
Lee ? ) 
- x = ne au 


yor at, : ’ q 2 - iy) 
Bite ek) ee a J 7 ee 
ey ee gilt’. | fp Riis hi ae. o 
reve Sd cg PIF 24 om ol — 
aah TC sa | 4 a 7 
% Po oa ee Py ; 
sart ead dee Ben ae eo ; ove > 
Viks Te mE + yee, aNTOD. 4 ee, gay 
. : h, 
* - a. et 5 


iit ea tk lac ‘% ‘ambos or f 
r. 3s bea p. * 1 - 


f ny ite : 
- sche eeetht jabs ‘pit oe Syantb ad a 


2 “ zi 


te ct ey lans. ft. oy bob a 4 aot 


FIGURE II-1 


CooRDINATE SYSTEMS USED IN ANALYSIS; SINK FLOW MopEL 
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The continuity equation for a symmetric flow field, in spherical 


coordinate system, becomes, 


r= (9) (11-1) 


ifs Ve = Vo = O (I1-2) 


where ¢ iS the cone semi-angle, Vo and Vo are the velocity components in 
spherical coordinate system. Normalizing the radial velocity OS) with 
respect to average velocity (V) in the capillary tube and radial position 
(r*) with respect to-the capillary tube diameter (D), and replacing f(o) 
by the constant Q* (32, 33) which is proportional to the volumetric flow 


rate (Q), Equation (II-1) can be written as: 


Vie ye pagal 
(“r*) ait (=) (11-3) 
V YD D 


Therefore, on logarithmic coordinates the slope of plot of the 
dimensionless radial velocity (Vix/V) versus the dimensionless radial 
position coordinate (r*/D), will be equal to -2 if Equation (II-3) is 
COYRect : 

In the study of Metzner, et al, Equation (II-3) applied to the 
central region where the cone semi-angle was less than or equal to 10°. 


In this region, Metzner, et al found that the velocity profiles at 


several axial locations upstream were nearly flat. 
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II.3  K-1 Model 

By taking photographs of the fluid column in Fano and fiber 
spinning flows and making diameter measurements, Kanel (26) found that 
the diameter of the thread was often given by: 
seal sae iy 74 aaa (11-4) 
dP 
in which Z is the axial position variable in a cylindrical coordinate 
system, d is the diameter of the fluid column, and p is a parameter. 


When p=1, it is the K-1] model, and Equation (II-4) can be written as: 


= ior 2. NC” (11-5) 


o.|— 


Considering the converging flow through a capillary tube, the 
axial position (Z) can be nermalized with respect to the diameter of the 


capillary tube (D). Then Equation (11-5) becomes: 


= Cie) 49 £2 (11-6) 


o.j— 


Figure (11-2) shows the coordinate system used for K-1 model in analysis 
of the velocity field of the converging flow into a capillary tube. 

From Equation (II-6) and the geometry in Figure (II-2), one can express 
the slope of the tangent to the outermost streamline at any axial 


location as; 


tan ¢ 5) 
C) 
2D 
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Figure 1-2 
CoorRDINATE SYSTEM USED IN ANALYsIS3 K-l MODEL 
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At the contraction Go 0), > becomes the initial cone semi-angle (¢ 0) 
and d becomes the diamter of the capillary tube (D). Then, 
C,D 

tan Pena Berea (11-8) 
Equation (II-7) shows that on logarithmic coordinates the slope of a 
plot of the tangent of the cone semi-angle versus the diameter of central 
core (d) will be equal to +2. 

Based on the flat velocity profile assumption by Kanel, the 


volumetric flow rate (Q) can be written as, 
2c 2 
21d evEG iS nd 
Q fo ae re ‘3 4 ) Vo 6¢ (TI-9) 


Substituting Equation (II-6) into Equation (11-9), and rearranging: 


We wblorl Ge SE 0g MeN aCze)é (11-10) 


Then, expanding Equation (11-10) in the form as: 
2 


2 2 
4 rc 
i dm teem ecy AF )Veleoi2e, Cate), <Aveati (11-11) 


Equation (II-11) shows a quadratic dependence of axial velocity on axial 
location (32,33), but this Equation is valid only for flat 
velocity profile. More formally, the volumetric flow rate (Q) can be 


expressed as: 
Garé=difs tances th Wide (11-12) 


where r is the radial coordinate in cylindrical coordinate system as 
shown in Figure (II-2). If the velocity profile is not flat, the axial 
velocity ey is a function of the radial location (r) and the depend- 


ence (cf II-9) is no longer on a. Therefore Equations (11-10) and (II-11) 
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16 
are no longer quadratics. However it is still of interest to test 
whether a quadratic dependence on axial location gives a reasonable fit 


and an equation similar to (II-11) will be assumed such that: 


Paver re Sty 4 el aa alto ie ade ero (11-13) 


The parameters K,, Ky and K3 can be found from the least square fit of 
the centerline velocity data. 
Differentiation of Equation (II-13) with respect to axial 


coordinate (z) shows the linear relationship of the dimensionless axial 


position coordinate (-) and the stretch rate along the centerline as: 
a Se 2K3) (Z é 
Vee 7 (BS + FY @ (11-14) 
II.4 Exponential Decay Model 


The velocity field within a melt spinning was analyzed by 
Kase (27) by solving the Navier-Stokes equations with the following 
assumptions: (i) steady state; (ii) constant density; (iii) axial 
symmetry both in temperature and velocity field; (iv) negligible gravity 
force; (v) negligible inertia forces; (vi) a know temperature profile. 
With the above assumptions, the viscosity was assumed known and was given 


by the expressions: 
Ho = pe’? (14 cr*) (11-15) 


where z and r are distances in cylindrical coordinates. Substituting this 
viscosity into the Navier-Stokes equations, a solution of the axial velocity 


(v,) have the expression: 


ty ard 


‘ede 7 


eee ae 


+9 sh ae e228 2 tah 


re RNA) ae ote ee, CO 


e 46 [4G oF dss 


ae 
. ‘ a be 
i“ 4 ry bi > ot + 
=| P| Dy 7 34) ul 
re e4 
Ceerr rs 
‘7 Ser tcd 
e ~~: 
; 
: Re 
ageye wows Pay “3, a 
reves © peer a0 
" oe ie nia a 
* - 
at PWT gt Bi . 
: ~~ ~~ 
% 
‘ + 
ry 54 j 
pat Pass shit bd, et 
ey Sex 
' | Ax 
a) hee t= eee 
ne y : LUOt, -& 
} 
=e 2 i rs 
} ) SAB SBS ; 
=| 
» P ] 
3 4) fy 
x Ls 
% 
b ots 2 * ty Hee 
‘ % rr cnet Tae 
me | } Bs 
; lsh 


= 7 - ‘ , ” <“¥ 
rit al ata de satgnitqooo ts 


vi “4 : ay f 6 hx wt n \ % o a tule ‘he att c . 
Av ' 


LZ 


i] 
V = V sen e(livia, to +a, netae Tr. > sac) © CLI-16} 


where be is assumed independent of r. 

The results of Kase showed that the velocity profile across 
the thread in any melt spinning could be assumed to be flat. In view of 
the generally high fluid viscosities and the absence of significant 
shearing, due to the low viscosity of the surrounding air, it would seem 
to be a reasonable assumption. In other words, it is reasonable to say 
that the values of ao, a,, ag -.-. are very small compared with unity. 


Therefore, Equation (II-16) can be written as: 
V = th jel ML Sig (11-17) 


Although the four parameters Lee ae Ho and g in Equation (II-15) and 
(11-16) are functions of axial position (z), they can safely be assumed 
to be constants by Kase in his study (27). 

Considering the converging flow through a capillary tube, the 
axial position (z) can be normalized with respect to the diameter of the 


capillary tube (D). Then Equation (11-17) becomes: 
a (4) 
TO ery ease (11-18) 


Equation (II-18) shows that on semi-logarithmic coordinates, a plot of 
the axial velocity (V,) versus the dimensionless axial location & 


will be a straight line with slope equal to a and intercept at 2 = 0 


equal to Vo" 
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Differentiation of Equation (II-18) with respect to axial 
position (z) yields: 


td geo G) 
a nan, e (II-19) 


Equation (II-19) shows that the stretch rate is an exponential function 
of the dimensionless axial location SL 
These three models will be used to analyze the results in this 


study and will be discussed in Chapter IV. 
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CHAPTER III 
EXPERIMENTAL EQUIPMENT AND PROCEDURE 


Lita Experimental Apparatus 
III.1-1 Entrance Region 

A reservoir of 2.64 cm internal width and 25.46 cm length was 
used as the entrance region. The internal length of the reservoir 
assured that the velocity distribution was fully developed prior to the 
location of the pressure tap which was located 10.20 cm from the 
anodized aluminium end plate of the reservoir. The reservoir was con- 
structed of 2.54 cm thick transparent lucite piate in order to allow for 
fiow visualization. 

Shown in Figure (III-1) is the device for holding the capillary 
tubes in place at the end of the reservoir. This holding device was 
6 


capable of withstanding maximum pressure of approximately 4.0 x 10 


rd i : 
dynes/cm” in the reservoir. 


III.1-2 Capillary Tubes 

Two glass capillary tubes with I.D. of 0.20, and 0.30 cm + 
0.00076 cm were used in this study. The capillary tubes were manu- 
factured by Fischer-Porter. The length to diameter ratios (L/D) for 
the two capillary tubes were 54.96 and 55.04 for I.D. of 0.20, and 0.30 cm 


respectively. 


III.1-3 Flow System 
An Instron (Model TT-BM) and a stainless steel hydraulic 
cylinder of approximately 15.60 cm I.D. were used to provide a constant 


flow rate of the polymer solution. The calibration for the Instron, 
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given in Appendix A, revealed that the flow rates calculated from knowledge 
of the Instron speed and the volumetrically measured flow rates were within 
0.5% of each other. 

When the piston was in its lowest position (the cylinder was 
nearly empty), there was still a clearance or space between the piston and 
the bottom of the cylinder. Before the cylinder was filled with polymer 
solution, a vacuum pump was used to eliminate air which would be trapped 
just below the piston. The air and any residual polymer solution in the 
space was drawn out by a vacuum pump into a flask. The See solution 
was then sucked into the cylinder by moving the piston upward. When the 
piston was in its highest position (the cylinder was full of polymer 
solution), it would stop automatically. This procedure was repeated 


each time for refilling the cylinder. 
BLI.2 Experimental Fluid 


I]I.2-1 Experimental Fluid 

The experimental fluid was a 0.2% aqueous solution of Separan 
AP-30. Separan AP-30 is a partially hydrolized polyacrylamide with 
molecular weight between 2 and 3 x 10° and was supplied by the Dow 
Chemical Company. The rheological properties of the polymer solution 
are given in Appendix B. 

A small amount of sodium bicarbonate was added to the polymer 
solution in order to protect it from bacterial degradation. Also, the 
polymer solution was used only once to ensure that the properties did 


not change from run to run owing to mechanical degradation. 
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File2=2z Tracer Particles 

Aluminium metal dust, manufactured by J.T. Baker Chemical Co., 
Lot number 2347, was used as tracer. Approximately 1 c.c. of these 
irregularly shaped particles were added to 300 1b. (136 Kg.) of the 
polymer solution. To prevent the particles from adhering to each other, 
they were dispersed with a glass rod in a small quantity of water 


(approximately 1 c.c.) before mixing with the polymer solution. 
bids Photographic Techniques 


III.3-1 Optical System 
The main features of the optical system which consisted of 
a light source, a converging lens, a rotating chopping disc, a reservoir, 


and a photographic assembly are shown in Figure (III-2). 


a. Light Source 
The light source used in this study was a CHRISTIE xenon 


arc-lamp which operated from a three phase CHRISTIE silicon rectifier 
unit with a maximum output of 3300 watts D.C. at 100 amps. The set of 
converging mirrors in the lamphouse (Model BSF 50) were adjusted to 
direct the light beam with its maximum brightness and light uniformity 
through a narrow slit located on the wall of the reservoir. The pro- 
cedures of adjusting those mirrors are shown in “Introductions for 
XENOLITE BSF Series Xenon Lamphouse". An adjustable horizontal slit 
placed immediately in front of the lamphouse lens, resulted in a light 
beam which was less in thickness than the slits on the chopping disc but 
wider than the approximately 0.41 mm wide slit on the wall of the 


reservoir. 
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b. Chopping Disc 


A paper disc, with 20 equal radial slits and spokes was 
clamped between two transparent lucite discs of 44.6 cm diameter. This 
chopping disc was located perpendicular to the light beam with its center 
in the plane of the light beam. The disc was driven by a variable speed 
drive in order to allow the time interval (the time for the slit to cut 
the beam of light) to be adjusted to a desired value. It is evident 
that higher flow rates required higher rotational speed in order to 
obtain reasonable lengths of streaks. The rotational speed of the 
chopping disc was measured by a stroboscope Type 1531-A STROBOTAC. 

Table (III-1) lists the R.P.M. of the chopping disc used 


according to different capillary tubes and flow rates. 
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Cs Photographic Assembly 

A Miranda 35 mm camera, located perpendicular to the light 
beam, was secured to Setco Building Blocks. This allowed the camera to 
be moved in planes parallel and perpendicular to the light beam. Three 
extension tubes, A. B. C. were used with a 1:2.8/105 mm lens and the 
resulting real magnifications on the film were approximately 0.73, 2.18, 
2.73 respectively. Tri-X film with ASA speed of 400 was used. Approxi- 
mate exposure times and apertures which were used in this study for each 


different extension tube are listed in Table (III-2). 


TABLE (III-2) 


EXPOSURE TIMES AND APERTURES 


Extension Tube Exposure Time Aperture 


used (sec) 


T4 
fio 
Tock 


L 


The exposure times listed in Table (III-2) were approximately 
the maximum exposure times before background light would cause loss of 
contrast or overexposure of the film. Therefore, actual exposure times 
were sometimes shorter than the values listed in Table (III-2), and were 


dependent on how much time was necessary until a few bright particles 


passed through the field of view. 
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oF Reservoir 

To minimize reflections and aid in providing as much contrast 
as possible for photographing the streaks, the reservoir was painted black 
inside except for the light slits on the side and the top camera viewing 
area. Also, the portion of the capillary tube outside the reservoir was 
covered by putting a sheet of black paper in front of it. During an 
experimental run the room was in total darkness. 

The upstream end of the capillary tube was made to project 
slightly into the upstream reservoir, so that the detailed velocity field 
could be studied photographically very near the capillary tube entry 
without the walls of the reservoir obscuring the view. The length of 
the capillary tubes protruding into the reservoir from the end plate were 
approximately 0.10 and 0.20 cm for tube I.D. of 0.20 and 0.30 cm 


respectively. 


III.3-2 Streak Photographs 
| Over 1,000 photographs were taken for three different flow 

rates with each capillary tube. Typical streak photographs are shown 
in Figures (III-3), (III-4) and (III-5). 

The flow field shown in Figure (III-3), was obtained with the 
0.20 cm I.D. capillary tube at the largest flow rate of 0.6392 cm/sec 
(flow moving from bottom to top) for which the velocity profiles were 
measured. The film exposure time was 8 sec. and the extension tube A 
was used. The time interval for a streak was 0.00787 sec. This figure 
shows the typical wine glass shaped core surrounded by large eddies. 

The flow field shown in Figure (III-4), was obtained with the 


0.20 cm I.D. capillary tube at a lowest flow rate of 0.1598 cm/sec 
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FIGURE III-3 Streak Photograph, 0.20 cm I.D. 


Q = 0.6392 cm?/sec, Extension Tube A 
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FIGURE III-4 Streak Photograph, 0.20 cm I.D. 


Q = 0.1598 cm?/sec, Extension Tube C 
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FIGURE IJI-5 


Streak Photograph, 0.30 cm I.D. 


Q = 0.1598 cm?/sec, Extension Tube B 
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(flow moving from bottom to top). The film exposure time was 30 sec and 
the extension tube C was used. The time interval for a streak was 0.00987 
sec. The flow field shown in Figure (III-5) was obtained with the 0.30 cm 
I.D. capillary tube at a lowest flow rate of 0.1598 cm?/sec (flow moving 
from bottom to top). The film exposure time was 25 sec and the extension 
tube B was used. The time interval for a streak was 0.0429 sec. Those 
two figuresshow a better view of the velocity field very oe to the 
tube entry than Figure (III-3). 

Using several photographs similar to those illustrated by the 
figures, the streamlines and the axial velocity distributions were deter- 


mined for several flow rates. 


III .3-3 Analysis of Streak Photographs 
a. Local Axial Velocity Measurements 

A Bell & Howell projector (Model 745C) was used to project 
the photograph negatives on graph paper which was graduated in 0.10 
inch divisions. The streaks were traced directly onto the graph paper. 
For example, consider the data obtained for the 0.30 cm I.D. capillary 
tube and for a flow rate of 0.3196 cm?/sec. Figure (III-6) illustrates 
how streaks would look as traced onto the graph paper. Figure (III-6) 
is in scale with respect to the tube diameter except the streaks of 
streamline number 11, which has been included to illustrate a streak 
close to the outermost streamline. In this example, the actual magnifi- 
cation used on the graph paper was 49.8. Knowing the magnification and 
the time interval, these streak lengths were converted into velocities 
along the streamlines. Although not precisely the axial or "z" 
velocity the small cone angles resulted in only small difference between 


axial and measured velocity. Therefore, these measured velocities can 
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be reasonably assumed as the axial velocities all through this study. 


Then, 


ye & (111-1) 
Mt 
where AL = _ length of streak, cm 

le = axial velocity, cm/sec 

MM = magnification of the streaks on the graph paper 
t =  ( 60/Nw), sec 
N = number of slits and spokes on the chopping disc 
7) = number of revloutions per minute, R.P.M. 


A sample calculation will be shown in Appendix F. 

To minimize the error due to the haif-shadow at the beginning 
and the end of every light period, and also the error due to finite 
particle size, nearly all the streak lengths were measured from leading 
edge to leading edge of the particles. Near the entry and close to the 
centerline of the capillary tube where the streaks were the longest, this 
was not always possible. 

The resulting velocities as a function of distance from the 
entrance are plotted in Figure (III-7) for the streamlines corresponding 
to those in Figure (III-6). The solid lines in Figure (III-7) are the 
eye-fit lines to the data and which have been extrapolated to estimate 


velocity right at the contraction. 
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Ficure I1I-6 
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b. Cone Angle Measurements 


The cone angles at the contraction and upstream of the 
contraction, are defined as the angle formed by the line tangent to the 
outermost streamline that enters into the capillary tube. This stream- 
line however, was not precisely determinable for all of the experimental 
runs. For example, in Figure (III-6), it is difficult to tell whether 
the streamline number 11 is exactly the outermost streamline of the flow 
that enter into the capillary tube or not. By using streamline number 


11, however, the initial cone angle was estimated to be Gree. 


Cx Discussion of Error 
In determination of the velocity distribution there were 


several sources of error associated with the projection technique. 


(i) At speeds of the chopping disc below 110 R.P.M., the 
speeds could only be measured by means of multiple images 
using a stroboscope. Repeat readings gave deviations as 


large as + 2% from the average. 


(ii) Using the projection technique it is important to know 

how accurately the streaks could be located with respect 

to the centerline of the capillary tube. This is because 
for each velocity profile, the data points were measured 
from several separate photographs. There were two methods 
used to locate the centerline of the capillary tube on the 
graph paper depending on which extension tube was being used. 
For photographs taken by extension tubes B and C where the 


walls of the reservoir could not be seen, the centerline of 
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the capillary tube was determined by measuring half way 
between the two visible outside walls of the capillary tube. 
By knowing the outside diameter of the capillary tube, the 
magnification of the photograph project onto the graph paper 
could also be found. For photographs taken by extension 
tube A, where the walls of the reservoir could be seen, the 
centerline of the capillary tube was determined by finding 
the center between the two reservoir walls. The magnifi- 
cation in this case could also be found by knowing the out- 
side diamter of the capillary tube or, alternatively, from 
the known inside width of the reservoir. 

Although the capillary tube was located following the 
procedure of centering of capillary tube which will be 
discussed in Section III.4-2, perfect symmetry was not 
always achieved in the velocity profiles because of non- 
uniformity in the wall thickness of the tubes. This would 
also cause an error in estimation of magnification. For 
example, consider data for the 0.30 cm I.D. capillary tube 
at the highest fiow rate 0.6392 cn/sec for which the photo- 
graphs were taken by extension tubes A and B. The axial 
velocity profiles are shown in Ficure (IV-2), and it can be 
seen that the local velocities on the right of the center are 
somewhat higher than those on the left. Also the experi- 
mental flow rates which were calculated from the areas under 
the curves of plot Ve versus re at different axial positions, 


were always less than the calibrated flow rate by about 6%. 
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Due to hazy beginning and end of the streaks, there is an 
error in measurement of the length of the streaks. This 

is more significant for short streaks than the long streaks. 
The rotational speed of the chopping disc was chosen to have 
most accuracy for streaks near the entry of the capillary 
tube. Therefore, streaks at locations further upstream 
were sometimes too short to measure with high accuracy. A 
good example for this case is for the 0.20 cm 1.D. capillary 
tube at a flow rate of 0.6392 cm?/sec and at an axial 
location equal to three tube diameters (z/D = 3) which was 
the furthest upstream velocity profile measured. Here an 
error of 12% between the experimental flow rate and the 


calibrated flow rate was found. 


Since the depth of field of the camera was always greater 
than the depth of the beam of light, it was not possible to 
distinquish whether the aluminium dust particles were at the 
top of the beam of light (nearest to the camera), or at the 
bottom of the beam of light or at any intermediate plane 
Tying within these limits. Thus, the streamlines on the 
photographs were not necessarily in the plane of the center- 
line of the capillary tube. The larger the velocity 
gradient the larger the possible error. Assuming the 
aluminium dust particles are randomly distributed on the 


fluid, this would cause a random scatter in the data. 
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(v). Since the measured length of streaks were not actually the 
distance travelled in the z-direction, the local velocities 


2 ony 
reported are not V, but (V_ > + Vp 73, This error becomes 


Z 
more significant when the angle between the streak and the 
centerline of the capillary tube becomes larger. Considering 
the worst case in this study, for the 0.30 cm I.D. capillary 
tube at lowest flow rate and at axial location equal to two 
diameters (z/D = 2), the derivation of (v, + Ve) from Up 
was less than 10%. While at higher flow rate and at axial 
location near the entry of the capillary tube, this 
derivation was always negligible. Therefore, no distinction 


Z 


2 L 
is made between es and (V, + Vy. 2 in this study. 


Thr4 Procedure 


TII.4-1 Preparation of Polymer Solution 

To prepare the 0.2% (by weight) Separan AP-30 Solution, a 
stainless steel tank was filled with 300 lbs (136 Kg.) tap water and 
272.1 gm. of solute. The salute was added slowly to the water so as to 
ensure that the powder was wetted without formation of large lumps and 
entrained bubbles. At the same time, a mixer agitated the solution very 
gently to avoid mechanical degradation and continued to do so for a 
period of approximately 10 hours after the final amount of solute was 
added to the solution. The solution was then left unsheared overnight. 
The specific gravity of the polymer solution determined for a fresh 


solution at 23°C, was 0.9984. 
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III.4-2 Centering of Capillary Tube 

Since the wall thickness of each capillary tube was not 
uniform.it was necessary to center each capillary tube so that the center 
of the capillary tube and the center of the light beam were on the same 
horizontal plane. 

The empty reservoir was first adjusted so that it was in 
a horizontal plane. The capillary tube was then inserted into the 
holding device which was attached to the end plate of the reservoir. 
The light source was turned on and the room lights were turned off. 
The light beam was then clearly visible against the aluminium end wall 
of the reservoir and the side of the reservoir opposite the side where 
the adjustable slit was attached. By adjusting the level of the light 
source a plane passing through the center of the light beam was in line 
with a plane passing through an etched mark located on both sides of the 
reservoir. These etched marks represented a plane passing through the 
center of the hole in the aluminium end plate where the capillary tube 
was to be located. The capillary tube was rotated by hand until the 
center of the beam of light and the center of the capillary tube were 


judged by eye to be in the same plane. 


IIT .4-3 Focusing on the Center of the Capillary Tube 

To focus on the center of the capillary tube a wire 
sharpened to a point on the end and having a diameter approximately 
equal to the inside diameter of the capillary tube, was inserted into 
the tube until the end of it protruded into the reservoir. The light 
source was then turned on and the shadow of the wire was clearly visible 


against the side of the reservoir opposite the side where the adjustable 


_ sitet utes ese erin ond to wobpt 


dipti oft |. viowraest ‘ads ‘ont € 


B 225 
r 72 rs 


or 
o ps +e — > 7  & - a 
Tiovvses) 2dG° 76 sista 
= 7 
=o iocrvard “te ce sitoert 
; . 
one mune s ond .sen rege 


i 
= Tel Wy eh Ls oT ¢ a f af 240 "I iv 
re 


ao f 
a 


~ 4 
Osh 


teutbs a de srstin 


ae eee eo 
ale ao ese . ~ r; a oe 
bu i 26? G2 oS32 


4 


“ns 79 +8voi Ooms wn (teh Ee cas 
e + “9g + 


‘ie 


| ‘gy mbit muta. 4 pent | 

| | a ec 

iene baat vd bagstor anw ath ROT te 152 oat et ate 630) 
- 


iTtas> sft te ncn 99 ats dns ogi f % aed st: 
i aut 
‘ntlest mtg odd AP ed oh dys 


‘ a ni. 7 - i eae yf: 
adust 5) ott te eden) we yr 
aren - ; , s Ps al n 
siiw 6 atut yall igs ao ad a A 


S| ~~ 
KONE. 2 ramet A Loaharyt sid 


baz : sent eon 9g ‘ont 


<< 


a 


ne. bh ie y ve 
7 ake ay 7 La , bats 


eu 
Dan f 


0 ae ee 


4] 


slit was attached. The wire was rotated slowly by hand until the shadow 
of the tip of the wire on the side of the reservoir was judged by eye to 
be in the center plane of the light beam, in other words, in the center 
plane of tne capillary tube. The camera was then focused onto the 
illuminated tip of the wire. 

Since the depth of field of the camera was greater than the 
thickness of the light beam, all aluminium particles within the beam's 


thickness were illuminated and recorded on the film. 


III .4-4 Data Collection of Velocity Profiles 
Measurement of the velocity profiles for the 0.2% weight 
solution of Separan AP-30 was discussed in the previous sections of this 


chapter. The results are tabulated in Appendix C. 
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CHAPTER IV 
RESULTS AND DISCUSSIONS 


Since the unstable velocity field of the converging fiow of 
viscoelastic fluid is nearly impossible to measure, only the velocity 
field for stable conditions is measured in this study. The highest 
flow rate chosen in this study so as to measure the velocity profiles 
within the central cord under stable conditions is 0.6392 cm/sec for 
both capillary tubes. Although the critical flow rates for the break- 
down of the entrance velocity fields for different capillary tubes 
estimated by Catania (10) was approximately 0.318 cm?/sec fOGea O32. em 
capillary tube, the velocity field at 0.6392 cm/sec for O0-c cm 1 Ds 
capillary tube is observed to be stable in this study. In this study. 
the axial velocity profiles at different locations upstream of the 
contraction were measured for three different flow rates of 0.6392, 
0.3196, 0.1598 cm?/sec for both capillary tubes. In addition, the 
initial cone semi-angles and the upstream cone semi-angles were measured 
so as to present a complete description of the flow field within the 


central cord of the converging flow. 


IV. Axial Velocity Profiles 

Table (IV-1) shows the flow rates ane the dimensionless 
locations upstream of the contraction where velocity profiles were 
measured. As indicated by the table the velocity profiles were measured 
in the region very close to the contraction. This region is of 
particular interest because it has been realized that the distortion of 


emerging extrudates may be closely related to the flow in this region. 
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TABLE AV=1 


Flow Rates and Axial Locations for Velocity Profiles Measured 


—_— eS eS Oe ee eee See Oe ee 


ees Flow Rate Axial Location (z/D) 
(cm) (cm?/sec) (No. of Diameter) 


0.20 0.1598 0.40% 
1.40 

2.0 

0.3196 40% 

150 

2.0 


(an) w 
(2, 
% 


0.6392 


0.30 0.1598 iets 


0.3196 
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TABLE IV-1 (continued) 


LE CD I TS SR tc i ce i tA i et ttt tests ihm tansicsdhciacintndnemets 
a a a ee ee ee a ee ee LD RIE LIE TELE, 


iO; Flow Rate Axial Location (z/D) 
(cm) (cm?/sec) (No. of Dia.) 


0.30 0.6392 0.0* 
1.0 
20 
0 


* From eye fit extrapolated data 


Raw velocity data are tabulated in Appendix C. From these 
data, typical velocity profiles are given in Figures (IV-1) and (IV-2) 
at several upstream axial locations for two different capillary tubes. 
The accuracy of these data may be judged in the following way. The 
local velocities can be integrated over the cross section of the con- 
verging flow that enters the tube and the result compared with the known 
volumetric flow rate. This comparison is of interest as it can be made 
in the upstream field for different locations. The mean deviation of 
these integrated flow rates from calibrated flow rates was approximately 
equal to 5%. 

Figure (IV-1) shows the velocity profiles which were obtained 
with the 0.20 cm I.D. capillary tube at the smallest flow rate 0.1598 
cm?/sec. Figure (IV-2) shows the velocity profiles which were obtained 
with the 0.30 cm I.D. capillary tube at the highest flow rate 0.6392 
cm/sec. The axial velocity profiles at the contraction, z/D = 0.0, 


were impossible to measure directly in this study but were obtained with 
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Ficure IV-1 


AxiaL Vevocity Prorives, 0.20 c4 1.D, Capittary Tupe, Q = 0.1593 om?/sEC 
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Ficure IV-1 (continueD) 0.20 co I.D. Caprttary TuBe 
Q = 0.15% ar/sec 


Vz (cm/sec) 


Vz (om/sec) 
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Ficure IV-2 
Axia Vetocity Prorites, 0.30 cm I.D. Capitcary Tupe, Q = 0.6392 c11/sec 
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Ficure IV-2 (CONTINUED) 0.30 cm I.D. Caprttary Ture 
Q = 0,6392 a?/sec 
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the extrapolated data points. This extrapolation procedure will be 
discussed in Appendix C. 

All of the data as illustrated in Figures (IV-1) and (IV-2) 
show that the velocity profiles within the central core are not flat as 
assumed by many investigators (26, 27, 32). In Figure (IV-2), for 
example, at z/D = 1.0, the velocity at the centerline and the velocity 
at outermost streamline were 10.9 cm/sec and 2.30 cm/sec respectively. 
Also, for the entire range of the flow rates in this study, the ratio of 
the centerline velocity at the contraction to the centerline velocity for 
developed flow of a power-law fluid in the capillary tube, was approximately 
constant and equal to 0.946. 

A plot of centerline velocities at the contraction (Vz,¢,0) 
which are tabulated in Appendix C versus the average velocity at the 
contraction (V) is shown in Figure (IV-3). The solid line with slope 


equal to 1.0 can be written as Equation (IV-1); 


ae AY (IV-1) 


With a power law fully developed velocity profile, the axial velocity of 


- the fully developed profile in the capillary tube can be expressed as: 


paces (n+1)/ 
eal. Ge) fee (1V-2) 


Then the centerline velocity of the fully developed profile is, 


3n+] 
n+] 


V fully devel. = ( ) V (IV-3) 


2a 


where n is the flow behavior index and equal to 0.55 from Appendix B. 


Substituting Equation (IV-3) into Equation (IV-1): 


ks = 0.946 \,, 


ee fully devel. (1V-4) 
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Ficure [V-3 
AxTAL CENTERLINE VELOCITY AT CONTRACTION VERSUS AVERAGE VELOCITY AT CONTRACTION 
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The ratio, Lo ey devel.) indicates the degree to which the 
velocity profile was developed at the contraction. Clearly the velocity 
profile is very nearly fully developed before the fluid enters the 
capillary tube. 

It is also of interest to compare the wall shear stress at 
entrance with those of fully developed tube flow. Considering a power- 


law fluid, the wall shear stress (x) can be expressed as, 


n 
ty = ko (S32) (1V-5) 


: rR 
Differentiating Equation (IV-2) with respect to r and substituting 


into Equation (IV-5), and rearranging, 


n n 
e 3n+] 8V. 
Tw or Ko ( An ) ( D? 
n 
= ko! (BY) (1V-6) 
‘tel 3n+1," ’ ; 
where Ko = Ko (aie - the values of Ko’ can be obtained from 


Equation (B-6) in Appendix B and equal to 2.24 (dynes/cm)-sec9*>”, 


Therefore the wall shear stress (11) of fully developed tube flow can be 
calculated from the following equation, 


0-55 


¥ 8V £ 
Aen 2a) (IV-7) 


The wall shear stress at entrance can be determined by Equation (IV-5) 
where (2) ep can be estimated from the slope of the tangent (at r=R) 
to the velocity profile at the contraction. The results of the wall 
shear stress at entrance and the wall shear stress of fully developed 


tube flow are tabulated in Table (IV-2) and are plotted in Figure (IV-4). 
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Ficure IV-4 Watt SHear Stress at ContRAcTion Versus 
Watt SHEAR Stress oF FuLty Devecopep Tupe FLow 
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The solid line in Figure (IV-4) can be represented as Equation (IV-8). 
(c) at entrance = 0.850 (t,) fully devel. (IV-8) 


The above equation also shows that the wall shear stress at 
entrance of the capillary is very close to the wall shear stress of the 


fully developed tube flow. 


TABLE IV-2 
Comparison of Wall Shear Stress at Contraction with Wall 


Shear Stress of Fully Developed Tube Flow 


r-D. Flow Rate (x Jat entrance | (r )fully devel. 

(cm) (cm/sec) (dynes/cm?) (dynes/cm*) 

0.20 0.1598 36.4] 41.69 
0.3196 Syeieae te} 61.0] 
0.6392 78.38 89.35 

0.30 0.1598 Th a7 21.34 
0.3196 20a ae 
0.6392 Shaoc 45.75 

IV.2 Cone Angles 


The procedures for measurements of the cone angles were dis- 
cussed in Chapter III. The measured cone semi-angles are tabulated in 


Appendix D. 
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IV.2-1] Initial Cone Semi-Angles 

By using the Oldroyd-Maxwell constitutive model and assuming 
a sink flow model of the velocity field within the central core of the 
converging flow to obtain the physical components of the normal stresses, 
Metzner, et al (32) and Murch (33) suggested a dimensionless group 
(r*°/0-, Q*) which is the reciprocal of the relaxation time - deformation 
rate product. The deformation rate (Q*/r*?) at the orifice was given 
by Metzner (31) as (Q*/R°) which was in terms of the orifice radius (R). 
By using the diameter of the capillary tube (D) instead of the radius (R), 
a plot of initial cone semi-angle ($0) versus a dimensionless group 
(D°/e,Q) on the logarithmic coordinates is shown in Figure (IV-5) 
60: the relaxation time at zero shear, is assumed equal to 1.0 sec (10). 
In Figure (IV-5), two parallel straight lines with slopes equal to 0.6 
are found for the 0.20 and 0.30 cm I.D. capillary tubes. The data in 
Figure (IV-5) can be reasonably represented by Equations (IV-9) and 


(IV-10) which are shown as the solid lines. 


ps, 0-6 
ee pera eee for 0.20 cm I.D. (IV-9) 
o 6 Q 
ps, 0-6 
6, = 27.20 (sq) fOrs@ecor eualey- (IV-10) 
0 


where $6 is in degrees. Considering Equations (IV-9) and (IV-10) and 
rearranging them, one can obtain an equation in the form of: 
D 0.6 
ee (Dy Gy) (IV-11) 
where V is the average velocity at the contraction and f(D) is a 


function of the capillary tube diameter (D). 
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Figure IV-5 Initia Cone Semi-AncLe Versus ae 
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Oliver and Bragg (38) suggested that a single equation could 
not be expected to give exact cone semi-angles for a range of visco- 
elastic solutions. However, for one solution and a limited range of 
Operating variables this might be possible. An empirical equation for 


initial cone semi-angle (¢,) of the form: 


i] 6 
Sino mom WK oey (IV-12) 
was suggested by Oliver and Bragg. Considering the small initial cone 
semi-angles in Oliver's study for low flow rates (¢, < Tee where the 
initial cone semi-angles increase when the flow rate decreases, it is 


reasonable to rewrite Equation (IV-12) as: 


Va aia (IV-13) 
sin ¢ 
% 


since the ratio ( ~\ approaches unity when ¢, is a small angle. 


Comparing Equations (IV-11) and (IV-13), these two equations 
are the same with f(D) = K ae (where t = -0.6). Therefore, a plot 
of initial cone semi-angles ($0) versus the average velocity at the 
contraction (V) should form a set of straight lines with the same slope 
equal to -0.6. The plots are shown in Figure (IV-6). Equations 
(IV-14) and (IV-15) are shown as the solid lines in Figure (IV-6). 
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Mie= .19s72 Vio” efor Os20cemar sD: (IV-14) 
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15.31 V" for 0.30 cm I.D. (IV-15) 
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Ficure IV-6 Initrat Cone Semi-Ancue Versus AVERAGE VELOCITY AT CONTRACTION 


100r 
Symbol I.D. (cm) 
(o) 0.30 
@ 0.20 
10 
re 
G 
@ 
ae -0.6 
mm ey AY 
rab) 1@) 
5, 
ee” 
1) 
~~ 
<r, = 9, 72v 08 


lOO =—200 


0:1 : 
iO 
V (cm/sec) 


WOTTaAaTyou 


prmwteem pent Soa ‘a ett daallaticnes ts crates 
7 | 
‘ 
} - 
¥ Be 
; { chi CF, T ; ke 4-3 
Soe - hat bo 
eo 
4 *4 ‘+ ar 
. ¥ 
G i. c 
{ v f 
: 
: 
i we zi A ba 
5 
es . ts i 2 
. be = 
rl z 
al co 
t . a 
t 7 
t ae , § 
y ok a 
i Fi ; ne 
: ee Py Ie me ve ad i . re 
: a 
4 ee 
Vi£.ci » 5 
v 
4 * -. af fl 
‘ od 
: t 
a ; ~ ne 
‘ emia eee) « 5 “Ka 7 
om ? reek * @ a & au 
: a s 
= i = 
; i - 
H Pa 
, - : 
: 
' le - 
a ; 
' 
} 
; . — na _ 
} ' 
' PN 
‘ “ 
% ' 
4 « . 7 
: “s va 
; : , | ) 2 
ER Ae Me) ay Se A Cotes Meena " 
Pr <> farye Pa 6 i” Ly : at pre 
wis wis. ba? t bp o : , Sie 
P Teoh ' 
s ~ ie, 
a\, “y ¥ re lh 
») Beene 
4 ' : Si Mee! a) Veal 
ey et ee 7 oe ae \¥' ( 


ig] aa ba ee ni 
i: Ce a ae 
hn f wi AT 


58 
The values of K and s in Equation (IV-13) can be easily found 
from Equations (IV-9) and (IV-10) or from Equations (IV-14) and (IV-15). 


Then an empirical equation: 


WRC Ree wie 


oO 


=0.0 


V (IV-16) 


can approximately represent the initial cone semi-angle for the 0.2% 
Separan AP-30 solution. 


From Equation (IV-16) the group (pl 1 ,y0-6) 


or approximately 
(D2/V) is suggested to represent the relationships between the velocity 
field within the control core, and the capillary tube diameter and flow 
rate instead of the dimensionless group (D°/e,Q). Figure (IV-7) shows 
a plot of the initial cone semi-angie (¢. ) versus (D°/V) on logarithmic 


coordinates. The data in Figure (IV-7) can be representea by Equation 


(IV-17). 
pe, 0-6 | 
¢ = 64.05 (7-) (IV-17) 


Thus for 0.2% Separan AP-30, the initial cone semi-angle is 
approximately represented by Equation (IV-17) for the values of Oy 
between 0.0020 and 0.0398 cm sec. The initial cone semi-angles cal- 
culated from Equation (IV-17) have an average error less than 3% of the 
measured values. 

The initial cone semi-angles measured by Oliver (38) for 0.1% 
ET 597 solution at low flow rate, where the initial cone semi-angles 
increase when the flow rate decrease, are also plotted in Figure (IV-7). 
Those Oliver's data can be reasonably represented by ¢ = 13.95 (pe/v)or4. 


The values of s and t in Equation (IV-13) corresponding to the Oliver's 


condition are 0.8 and -0.4 respectively which are not the same as the 
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Ficure IV-7 Initia. Cone Semi-ANcLE Versus (D2). 
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values in this study. This result indicates that the values of s and t 
in Equation (IV-13) are constant only for one concentration of a polymer 
solution in the limit range of operating variables, they will change for 
different concentrations, different polymer solutions, or different range 
of operating variables as mentioned by Oliver (38). 

Unfortunately, the group (D°/V) in Equation (IV-17) is not 
dimensionless. But being an empirical equation determined by the experi- 
ment data in this study, Equation (IV-17) is still valid to represent the 
initial cone semi-angle in the limited range of operating variables in 
this study. A possible dimensionless group of these variables would be 
(o"7e “ve ) which is the ratio of Reynold number ee to the cube of 
Deborah number (e). From the data in this study, the value of this 
group was estimated approximately between 4x10°° and 4x10°°. ‘In view of 
the large radial dependence of velocity, the necessity of including shear- 


ing viscosity as an important variable is not surprising. 


lye2.2 Upstream Cone Semi-Angles Development 

Since the initial cone semi-angle (90) is a function of (D°/V), 
the cone semi-angle (4) are plotted in Figure (IV-8) to (IV-10) as a 
function of the same group so as to compare the development of the cone 
semi-angles for different I.D. of the capillary tube at the same flow 
rate. From the data plotted in Figures (IV-8) to (IV-10), several 
observations can be made as follows: 

(i) For some capillary tube, the gradient of the cone 

semi-angles along the axial locations increases as the 


flow rate decreases. 
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(ii) For the same flow rate, the gradient of the cone 
semi-angles along the axial locations increases 


as the diameter of the capillary tube increases. 


The upstream cone semi-angles development predicted by K-] 
model was shown in Chapter II, to be 
Cy 
ete dé 


tang = - (IV-18) 


2D 

On logarithmic coordinates the slope of a plot of the tangent of the 
cone semi-angle (tan o) versus the diameter of the central core (d) will 
be equal to +2 if Equation (IV-18) is correct. Figures (IV-11) and 
(IV-12) show the data of this study. The solid lines representing 
Equation (IV-18) has good agreement with the data only at axial 
locations further upstream from the entrance and give less than experi- 
mental values of the cone semi-angles near the contraction. Because 
of this poor agreement with the data near the contraction, the K-1 model 
is not suitable in representing the upstream cone semi-angles develop- 
ment in this study. 

The cone semi-angles and the diameters of the central core at 
several upstream axial locations are tabulated in Appendix D and serve 


to give a description of the shape of the central core. 


Nc Central Core Velocities 
To obtain an empirical description of the velocity profiles 
within the central core which would describe the dependence on axial 


location, the three models described in Chapter II were considered. 
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Ficure IV-l1 
UpsTREAM Cone Semi-Anctes Devecopment; K-1 Mopec; 0.20 cm I.D. CaPrLearY Tupe 


18) 

5 Q (cm?/sec) z/D for points shown 

0.1598 0.0) 1.03 ;2.0573.09 74.0, 5.0, 

0.3196 O.0 91. F..0 5.240 Fa sa0 AO fo, 

° 0.6392 0.09). 5..0 4.210 2e3'50 64,0 2540's 
| 
OS 
ce 
“O01 
0-05 
0-02 
0-0! 

OOl 0:02 OO5 (Ol Ec 0-5 | 2 


d (cm) 


neon # 


- 
0.8 0.2 04 
: 


va SE .G.e 0.2 0.1 
; 5 ' 

| 

£ 


5 a 
ft 
: 

t x ; 
' : - : 
© » 
, 1 
a zat ¥ 
<3 7 
of + 
a » 
Th 
z a 
for % - 
es = 
x f ’ 
i 
, 
Pre 
BR : 

% 
bs 
i 


a 
, ad + 
: 


fgg oto 


epee, £35 


at Pe "een 
rwar i: pie A: 
ae ot 


TEs 


2 
=) 
ean) 
zy 
iy 


oP 
cy 


Seat alll 
% 


66 


Ficure ]V-12 
UpsTREAM Cone Semi-ANcLe Devetopment; K-] MopeL; 0.20 cm I.D, CaPrLLarY Tupe 
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Only the centerline velocities were considered to obtain the relation- 
ship of the dependence on axial location. The centerline velocities at 


different axial locations for several flow rates are tabulated in 


Appendix C. 


IV.3-1 Sink Flow Model 

As described in Chapter II and shown by Figure (II-1), the 
velocity field within the converging region may be analyzed employing a 
spherical polar coordinate system. The origin is located within the tube 
at a point determined by extending the converging streamlines to their 
intersection. The location of the origins of the sink flow model as 


determined by the streamlines are tabulated in Table (IV-3). 


TABLE IV-3 
Location of Origin for Sink Flow Model 


130: Flow Rate Location of Origin (-z) 
(cm) (cm?/sec) (cm) 
0.20 0.1598 iLo5 
0.3196 2.49 
0.6392 3.56 
0.30 0.1598 0.9] 
0.3196 1238 
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By knowing the location of the origin of the sink, one 
can easily convert the cylindrical axial coordinate (z) into spherical 
radial coordinate (r*) because only spherical coordinates are used to 
analyze the sink flow model. 

Consider the centerline velocities, Equation (II-3) can be 


rewritten as: 


oes ee (e) (IV-19) 


Therefore, the slope of a plot of (Vz,¢/V) versus (r*/D) on the logarith- 
mic coordinates will be equal to -2. Figures (IV-13) and (IV-14) show 
the data of this study where the solid lines are the straight lines with 
Slopes equal to -2. The results have good agreement only near the 
entrance. This was also suggested by Kanel (26) and Murch (33). At 
larger upstream axial locations, the sink flow model predicts velocities 
that are higher than the experimental data. This is because at larger 
upstream axial locations the streamlines are curved and are not directed 
to the origin. This deviation increases when either the flow rate 
decreases or the diameter of the capillary tube increases. 

For 0.30 cm capillary tube, at a flow rate of 0.1598 cm/sec good 
agreement was for upstream axial locations as large as 1.0 diameter, 
while at 0.6392 cm?/sec agreement occurred for upstream axial locations 
as large as 1.5 diameters. For the 0.20 cm capillary tube, at a flow 
rate of 0.1598 cm?/sec, agreement occurred for upstream axial locations 
as large as 1.5 diameters, while at 0.6392 cm?/sec agreement was up to 


2.0 diameters. 
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Ficure IV-13 
CENTERLINE VeLocity; SINK Fiow MopeL; 0.20 cm J.D. Caprtiary TuBe 
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Ficure ]V-14 
CENTERLINE VeLocity; SINK FLow Move; 0.30 a4 1.D. CaPittary TuBe 
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Because of the poor quantitative as well as qualitative agree- 
ment with experiment data, the sink flow model is inadequate to represent 


the relationship between velocity along the centerline and axial location. 


IV.3-2 K-] Model 
As described in Chapter II, the K-1 model suggests a quadratic 


dependence of axial centerline velocity on upstream axial location as, 


Ee Sa ae (2) (IV-20) 


The three parameters K; , K» and K3 were found from the least square fit 
of the experiment data for different flow rates and diameter of capillary 
tubes. The values of these parameters are tabulated in Table (IV-4). 
From these values of K; , Ky and K3 , it is clear that the relationship 
ko“ = 4K,K3 required by Equation (II-11) is not satisfied. This 


demonstrates the inadequacy of the K-1 model. 


TABLE IV-4 


Parameters of K-1 Model 


Eau. Flow Rate Ky K> K3 
(cm) (cm?/sec) (cm/sec) (cm/sec) (cm/sec) 
0.20 0.1598 8.3760 -2.194] 0.1957 
0.3196 1625370 -2./086 0.1450 
0.6392 Bb2.3205 -4.3701 0.2009 
0.30 0.1598 3.3450 -1.7700 0.2600 
0.3196 7.0720 -2.9519 0.3671 


0.6392 14.5700 -4.1914 0.4286 
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Typical experimental centerline velocities are shown in Figure (IV-15) 
and compared with the least square lines. Considering the differentiation 
of Equation (IV-20) with respect to axial location (z), yields: 


K 
V of (ee 


we (IV-21) 


nN 
U 


Typical stretch rates along centerline which are tabulated in Appendix E 
are shown in Figure (IV-16) and compared with Equation (IV-21). 

Although the quadratic equation has good agreement with 
the centerline velocity data as shown in Figure (IV-15) and also with 
the stretch rates as shown in Figure (IV-16), it is not obvious how 
to obtain the relationships for the parameters Ky and K3 with flow rate 
and diameter of capillary tube (the ratio of K, to the average velocity 
at the contraction is approximately constant). In other words, when 
the flow rate or the diameter of the capillary tube is changed, one can- 
not predict what the changes of Ky and K3 will be. Therefore the K-1 


and/or quadratic equations are considered to be of limited interest. 


IV.3-3 Exponential Decay Model 
As mentioned in Chapter II, the velocities within the central 
core can be expressed as a exponential function of the axial locations. 


Consider the centerline velocity, Equation (II-18) can be rewritten as: 


a Z a 
is = Vo exp, (a (5) ) (IV-22) 


Equation (IV-22) predicts that on semi-logarithmic coordinates, a plot 
of the centerline axial velocity (V, ¢) versus the dimensionless axial 
location (z/D) will be a straight line with slope equal to a and inter- 


cept at z/D = 0 equal to Vo. 
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Figure IV-15 Typical CENTERLINE AxraAL VeLocity; K-] MopeL 
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Ficure IV-16 Typical CENTERLINE StretcH Rates K-1 Mopet 
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Figures (IV-17) and (IV-18) show the experimental data of this 
study with the solid lines determined by the least square method. Good 
agreement with all of the centerline velocity data is obvious. In 
addition, some data from the other streamlines where experiment data can 
be measured as close as 0.25 diameters for the 0.30 cm 1.D. capillary 
tube are plotted in Figure (IV-19) with solid lines again determined by 
the least square method. From Figures (IV-17) to (IV-19) the exponen- 
tial decay model has good agreement with all the data not only for the 
centerline velocities but also for those of the other streamlines. The 
experimental data measured by Metzner, et al (32) for a 0.5% Separan 
AP-30 solution flowing from an 18" reservoir into a 3.76 cm I.D. tube 
are also plotted on semi-logarithmic coordinates in Figure (IV-20). 
These data which were chosen from the two streamlines very close to the 
centerline (6 = 0°20’) for two different flow rates have good agreement 
with the exponential decay model too. 

The two parameters Vo and o in Equation (IV-22) which were 


determined by the least square method are tabulated in Table (IV-5). 
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Ficure IV-19 Typicat AXIAL VeLocity; ExpoNenTIAL Decay MopeL 


Sym bol 1. Deen). (cm/sec) 
V 0.30 0.6392 


A 0.30 0.3196 
oO 0.30 0.1598 
0 0.30 0.1598 
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TABLE IV~5 
Parameters of Exponential Decay Model 


ae=a==eeeeeeeelllelleellll lS 


ag He Flow Rate Vo o 

(cm) (cm?/sec) (cm/sec) 

0.20 0.1598 8.5071 -0.2902 
0.3196 l6251z6 -0.1826 
0.6392 32.6827 -0.1486 

0-30 0.1598 eee WAT) -0.7289 
0.3196 7.6916 -0.4814 
0.6392 14.6698 -0.3078 


A plot of the Vo in Table (IV-5) versus the average velocity 
at the contraction (V) is shown in Figure (IV-21), the solid line is 


represented by Equation (IV-23): 
Vo = 1.645 V (1V-23) 


Compare Equation (IV-23) with Equation (IV-1), the difference between 
the extrapolated values of the centerline velocities at the contraction 
from the exponential decay model and the eye fit method is less than 2%. 

Since, in Section (IV-2.1), the initial cone semi-angle can be 
represented by a group (-), this group was also used when trying to 
obtain the relationships between the parameter a and the flow rate and 
diameter of capillary tube. Figure (IV-22) shows a plot of the para- 
meter a versus (7) on logarithmic coordinates. The data in Figure 
(IV-22) can be reasonably represented by Equation (IV-24) which is shown 
as the solid line in Figure (IV-22): 


ne. 0.545 
a = = 4.04 (F-) (IV-24) 
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Ficure IV-21 
ParAMETER Vo oF Exponentiat Decay fopeL versus AVERAGE VeLocity aT CONTRACTION 


60 
Vo 
(cm/sEc) Symbol I.D. (cm) 

©) 0.30 
50 ® 0.20 
40 
«@——Vo = 1.6448V 

30 

20 

10 P 

O 10 20 30 40 


V (avsec) 


, waitae at m0) Tv! hy OAH 
¥ phon “ i teen ener ra on 


ee 


& 
a 


he 
- 


a ee ee ee 
J 


i 
. 
~ 5 
i h H 
+4 
ee - 
+ 
3 
aS 
i. 
¥. 
Rae 
; - 
el 
. 
i + 
= 
6 
= 
- 
i € 
ee ee 
? 
q 
- 
] 


Bz | " oe! 
: A a te mn iarae iene ee 


hay ae + wry. te os 


a 


fas 
Wi 


i's 


100 


Ficure IV-22 PARAMETER Of OF THE ExponeNTIAL Decay MopeL 


Symbol I.D. (cm) 


p2_ 0-545 
gee GO eas 05() 
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The mean deviation of the yalue of a calculated from Equation (IV-21) 
from the values of o in Table (IV-5) was less than 5%. 

Substituting Equations (IV-23) and (IV-24) into Equation (IV-22), 
and empirical equation in the form: 


pe 0.545 


V = 1.645V exp [-4.04 (7) (a7 (IV-25) 


256 


was found to represent the relationships for the centerline velocity with 
flow rate and capillary tube diameter. Figure (IV-23) shows the 
comparisons of Equation (IV-25) and the experimental data for two typical 
runs. The predicted values of Equation (IV-25) have a maximum deviation 
of about 5% and a mean deviation less than 3% from the experimental data. 
Differentiation of Equation (IV-25) with respect to axial 


location (z) yields: 


V ne 0.545 pe 0.545 4 
V2 z5¢ = - 6.646 (5) (7) exp [-4.04(7-) (5) J 
(IV-26) 


The comparisons of Equation (IV-26) and the experimental data are shown 

in Figure (IV-24) for two typical runs. The predicted yalues from 
Equation (IV-26) have a maximum deviation of about 9% and a mean deviation 
of about 3% from the experimental values. Also it is of interest that 
Equation (IV-26) predicts a maximum value for the stretch rate of the 
contraction. Since the stretch rate of the fully developed tube flow 
must be equal to zero, therefore the stretch rate must decrease and go 

to zero somewhere between the contraction and the location where the 


velocity profile is fully developed. 
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Ficure IV-23 
ComMPARISON OF CENTERLINE AXIAL VeLocITy with ExpoNENTIAL DECAY MobEL 
Vo 3 
eae Symbol I.D. (cm) Q(cm" /sec) 
SEC 
@) 0.30 0.3196 
lO e@ 0.20 0.1598 
“* 
> Vv -0.289(2) 
5 Va z,¢ = 8.373e D 
° 
Q 
, 
rey 
) 
Zz 
Vee = Pisses oe 
O 


La 


- 
anes . 2 
=! ~~ - 7 * _ 


a 


- 


+7 


4 
a 
- laa 
wa 
i. a 
E oat g) 02 =F 
¥ = 


: 


Cee. ce er oes 
. 


be 


7’ 

¥ 

: : 

‘ ' 

a aT aan | 

a : . Ah 

‘ ‘ ee f: . a 

ena wt oe mm ps 

x ac ; 

Sao, oe 
= > —_ 

M : 


ae 


ey Phares) oar 


ai 


Ficure [V-24 
ComPARISON OF CENTERLINE STRETCH PATE WITH EXPONENTIAL Decay MopEeL 
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The above experimental results reveal] that the velocity field 
within the central core can not be predicted simply by sink flow model 
or K-1] model. However, the exponential decay model has good agreement 
with the experimental centerline velocities. Besides, the experiment 
data show that the ratio of the square of the capillary tube diameter to 
the average velocity at the contraction (D°/y) can be used to reasonably 
represent the relationships between the velocity field, and the capillary 


tube diameter and flow rate. 


AL pr mt aha ~ : aes of vite peg 
a ] ; 7 <p cf t) an . ; v 
ae i 0 he neg eich ye Me 


ie te) a a aN 


. as Ry P 7 ; . “i : "ag 4 Talks 
2 a ee } ; n 
ted ' ae at ) ms Hi hs pars =e ; 
ny ‘ ieee 


ae 
: fir eae 


‘ 


iu ba i 7 
was itt teen: 
ie ney) 4 i) oi 
Bai) i 

beate a 
ae 


87 
CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 


V.] Conclusions 
The conclusions of this study can be summarized in the following 
manner: 
(i) The experimental velocity profiles for 0.2% aqueous solution 
Of Separan AP-30 within the central core showed considerable 
development upstream of the entrance. At the entrance, 
centerline velocities were approximately 96% of the centerline 
velocities of the fully developed tube fiow and the wall shear 
stresses were approximately 85% of the wall shear stresses of 


the fully developed tube flow. 


(ii) The sink flow model and the K-1 model failed to predict the 
centerline velocities within the central core because the sink 
flow model had poor agreement of the experiment data at large 
upstream locations from the contraction, Parameters in a 
quadratic equation suggested by K-1 model did not have obvious 


relationships with the flow rate and the capillary tube diameter. 


(iii) The initial cone semi-angle could be expressed as: 
2 0.6 
, DL 
Pty 64.05 (y ) (V-1) 
2 


where ¢, is in degrees, for values of (7) between 0.0020 and 
0.0398 cm sec. The initial cone semi-angles predicted by 
Equation (V-1) had a mean deviation of about 3% from the 


measured values. 
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(iv) The centerline axial velocities within the central core could 
be expressed as: 


- 0s545 


Vg = 1-646 V exp [-4.04(7 (4) (V-2) 


within the range of operating variables in this study. The 
centerline axial velocities predicted by Equation (V-2) had a 
maximum deviation of about 5% and a mean deviation less than 


3% from the experimental data. 


(v) The stretch rates along the centerline obtained by differen- 
tiating Equation (V-2) with respect to axial location (z) are 
given by: 


2709545 _ non 


) exp [-4.007) (2) 


V, /D 
Y( ; 


(V-3) 
The stretch rates predicted by Equation (V-3) had a maximum 
deviation of about 9% and a mean deviation of about 3% from the 


experimental values. 


V.2 Recommendations 
(i) Data in this study are based on only two capillary tube diameters, 
therefore additional data should be generated to test the 


validity of the diamter dependence implied by the correlations. 


(ii) In this study the empirical equation developed, predicts 
the stretch rate is a maximum at the contraction. On the 
otherhand, the stretch rate at a downstream location in the 
capillary tube where the velocity profile is fully developed 


must be equal to zero. This means the stretch rate must 
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decrease and go to zero somewhere between the contraction 
and location where the velocity profile is fully developed. 


More detailed study within this region is recommended. 


Since all of the empirical relationships in this study are 
restricted as to the range of operating variables, some 
experiments with other polymer solutions are recommended in 
order to compare with the predictions such as Equations (V-1) 


and (V-2) in this study. 
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NOMENCLATURE 
@o, ays Ags parameters defined by Equation (II-16) 
C parameters defined by Equation (11-15) 
Gis? 05 parameters defined by Equation (II-6) 
Casio. parameters defined by Equation (II-4) 
d diameter of the fluid column within the center core,cm 
D I.D. of capillary tube, cm 
Dp diameter of the plate, cm 
kK> torsion bar constant, dyne cm/micron 
K parameter defined by Equation (IV-13) 
K' parameter defined by Equation (IV-12) 
Ki, Ko, Kg K-1 parameters, Equation (II-13), cm/sec 
M magnification of the streaks on the graph paper 
n flow behavior index in power-law 
N number of slits and spokes on the chopping disc 
p parameter defined by Equation (II-4) 
Q volumetric flow rate, cm/sec 
o* volumetric flow rate function defined by Equation (II-3) 
r radial coordinate in cylindrical coordinate system, cm 
o radial coordinate in spherical coordinate system, cm 
R radius of capillary tube, cm 
Rp radius of the plate, cm 
S parameter defined by Equation (IV-13) 
s parameter defined by Equation (IV-12) 


S$; sensitivity of the torque transducer, micron /volt 
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parameter defined by Equation (IV-13) 

parameter defined by Equation (IV-12) 

time interval for streak, sec 

torgue, dyne cm 

average velocity at the contraction, cm/sec 

velocity at contraction defined by Equation (II-16),cm/sec 
radial velocity in cylindrical coordinate system, cm/sec 
velocity components in spherical coordinate system, cm/sec 
axial velocity in cylindrical coordinate system, cm/sec 
axial centerline velocity, cm/sec 

axial centerline velocity at the contraction, cm/sec 

axial stretch rate, sec”! 

axial stretch rate along centerline, sec’ | 


axial distance upstream of contraction, cm 


parameter defined by Equation (II-18) 
parameter defined by Equation (II-16) 
parameter defined by Equation (II-15) 
movement of torsion head transducer, volt 
relaxation time, sec 

zero shear relaxation time, sec 

viscosity, poise 

parameter defined by Equation (II-15),poise 
kinematic viscosity, cm/sec 


shear rate, sec 
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shear stress, dyne/ cm? 


wall shear stress, dyne/cm@ 

cone semi-angle, degrees 

initial cone semi-angle, degrees 

angle of viscometer cone, degrees 

number of revolutions of the chopping disc/minute, RPM 


angular rotation of the plate, rad./sec 
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APPENDIX A 
CALIBRATION OF INSTRON 


An Instron (Model TT-BM) and stainless steel hydraulic 
cylinder of approximately 15.60 cm I.D. were used to provide a constant 
flow rate of the polymer solution in this study. At Push-button Speed 
Selector settings of 0.2, 0.1 and 0.05 cm/min, the measured flow rates 
which were determined by the volume of the polymer solution collected 
over a known period of time were 0.6397, 0.3202 and 0.1603 an fsec 
respectively. Figure (A-1) shows a plot of the measured volumetric 
flow rate and the Instron head speed or speed of piston. The solid 
line which is an eyefit through the data results in a calibration factor 
for flow rate of 3.196 (cm?/sec)/(cm/min). The comparison of the cali- 
brated flow rate and the flow rate calculated from knowledge of the 


Instron head speed and I.D. of the cylinder is tabulated in Table A-1. 


TABLE A-1 


CALIBRATION OF FLOW RATE 


Instron Head Speed Calibration Flow Rate | Calculated Flow Rate 
(cm/min) (cm?/sec) (cm? /sec) 
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From Table (A-1), the flow rates calculated from knowledge 
of the Instron speed and calibrated flow rates were within 0.5% of each 


other. 
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APPENDIX B 


RHEOLOGICAL PROPERTIES OF THE POLYMER SOLUTION 


A WEISSENBERG Rheogoniometer was used to measure viscosity of 
the polymer solution. This Rheogoniometer was calibrated with two 
Standard oils of different viscosity. The theory of this Rheogonio- 
meter which was shown in the Instruction Manual of the Weissenberg 
Rheogoniometer (55) is also discussed here. 

The Weissenberg Rheogoniometer is basically a simple cone and 


plate viscometer. The shear stress is given by: 


re Sly 
2m, 
(B-1) 
tht a 
Dp 
where t = Shear stress, dynes/cm= 
T = SyArks = torque ,dynes cm 
S} = torque transducer sensitivity, micron/volt 
Ay = movement of torsion head transducer, volts 
kt = torsion bar constant, dyne cm/micron movement of the 
transducer 
Ro = radius of the plate, cm 
DD = diameter of the plate, cm 


The shear rate is given by, 


5 ]|CUe)| angular velocity in radians/sec 
tan yp 


"9 vind " 
| a ss nae eats 
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i.e. for the small angles of yp used: 


oO = 180 2 
7 OY 
(B-2) 
bg 360 
yt 
where o = shear rate, sec’ 
yw = angle of cone, degrees 
t = 2n/2, sec./rev. 
2 = angular rotation of the plate, rad./sec. 
With t and o measured apparent viscosity is calculated by: 
<8 e 
Hoe a (B-3) 


For all runs in this study only one plate and one torsion bar 


was used. For this plate and torsion bar, the Rheogoniometer parameters 


are: 
DD = 7.50 cm 
y pee. 10) 
ky = 0.1066 x 10° dyne cm/micron 
Sy = 1.008 micron/volt 
Gap Set = 39 microns 


Therefore, Equations (B-1) and (B-2) become, 


T ae 0-973 Ar dyne/cm® (B-4) 


o = S228 SNe. (B-5) 
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By setting a rotational speed of the motor which drove the 
plate, the movement of torsion head transducer (Ar) was recorded by a 
recorder. Using this knowledge, t and Ay can be found. Then the 
values of shear stress and shear rate can be calculated by substituting 
t and Ay into Equations (B-4) and (B-5). 

Plots of shear stress and, viscosity versus shear rate are 
given in Figures (B-1) and (B-2). In Figure (B-1), good agreements of 
the data with the calibrated solid lines of shear stress versus shear 
rate with slope unity and the calibrated constant viscosity horizontal. 
dashed lines obtained for the standard oils, indicate the Newtonian 
behavior of these fluids. In Figure (B-2) the solid line of shear 
stress versus shear rate can be represented by Equation (B-6) within the 


S 


range of shear rate from 5.0 x 10° to 1.10 x 10 sec”! which is of 


interest im this study. 
Tee =) Lee eee (B-6) 


The dashed line of apparent viscosity versus shear rate can be represented 


by Equation (B-7) 


-0.45 
oO 


ie = 2.24 (B-7) 


Equations (B-6) and (B-7) indicate power law behavior of the polymer 
solution with flow behavior index (n = 0.55) in the range of shear rates 


studied. 
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APPENDIX C 
AXIAL VELOCITY PROFILES AND CENTERLINE VELOCITIES 


Cal Axial Velocity Profiles 

Contained in this section are the axial velocity profiles at 
several upstream axial locations for various flow rates. The method 
used to measure the axial velocities and associated error were discussed 
in Chapter III. In addition, included in the table are the axial 
velocity profiles at the contraction (z/D = 0). These were determined 
by drawing eye-fit lines through the arithmetic plots of measured axial 
velocities as a function of distance from the entrance and extrapolat- 
ing the curves to the entrance. 

Although the measured velocity was not previsely the axial 
velocity (discussed in Chapter III10, the small cone angles resulted in 
only small difference between measured and true axial velocity.  There- 
fore, these measured velocities are tabulated as the axial velocities in 
Table (C-1) to (C-6). The radial location in Table (C-1) to (C-6) are 
the actual distance from the centerline of the capillary tube as dis- 
tinquished from apparent distances on graph paper. The radial locations 
marked with an "*" means that these radial locations are on the opposite 
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TABLE C-1 (continued) 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 


(No. of Dia.) (cm) (cm/sec) 
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TABLE C-2 
Diameter: 0220 ceme LD. 
Flow Rate: 0.3196 cm? /sec 
Axial Location (z/D) Radial Location (‘) Axial Velocity (Vz) 

(No. of Dia.) (cm) (cm/sec) 
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Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 
(hoe0F Dia.) (cm) (cm/sec) 
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TABLE C-2 (continued) 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 


(No. of Dia. ) (cm) (cm/sec) 
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TABLE C-3 (continued) 
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Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
130 0.003% 26.16 
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TABLE C-3 (continued) 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 


(No. of Dia.) (cm) (cm/sec) 
2.0 0; 101+ 9.20 
0. WO 4.83 

3.0 0.000 21.10 
0. 01e* 19.60 

O2017= 19.60 

0.013 1oe50 

0.006* 19.25 

05025 18.80 

0.039% 18.70 

0.041] 18.50 

G.031 18.43 

0.040* 18.30 

0.046* WSS) 

0.064* 16.79 

0.059 16.54 

0.072* loc72 

0.066 15.60 

0.061* 15.60 

0.088 13.00 

0.083 13.00 

0.083* 12.53 

0.090 TE.50 


COnmiuedee ie 


by pee : + nt ad 


oC. 
i 
4, 8t 


ate 

cy.) Sepia 
bh .8r 
0c at 
ze.NT 
evar 
pavet . 
a =. 
ao.er 
03.05 
W.£f 
00. ef hs gen 

ofl 9. ae 
02 RF 


iv ‘ or 7 iy i mis 
stew a eee yi, ae 
7 2 7 a) | ie | 


Axial Location (z/D) 
(No. of Dia.) 


34.0 


TABLE C-3 (continued) 
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(cm/sec) 


-10 


125 


TABLE C-4 


Diametepre—0290-cmeIeD; 
Flow Rate: 0.1598 Gn feet 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 


(No. of Dia.) (cm) (cm/sec) 
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Radial Location (r) 
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Axial Velocity (Vz) 


(No. of Dia.) (cm) (cm/sec) 
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TABLE C-4 (continued) 
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Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
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TABLE C-4 (continued) 
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Axial Location (z/D Radial Location (r) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
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TABLE C-5 


Diameter: 0.307 Cmte Le 


Flow Rate: 0.3196 aie see 


—— 


Axial Location (z/D) Radial Location (Tr) Axial Velocity (Vz) 


(No. of Dia.) (cm) (cm/sec) 
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Radial Location (r) 
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Axial Velocity (Vz) 
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TABLE C-5 (continued) 


Axial Location (z/D) Radial Location (r) ~ Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
1.0 0.146% 2.30 
051.58" 2.20 
0.176% Ba70 
0.186 — 1.29 
220 0.000 2.97 
0.017 2.94 
0.009* Va Bas 
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0.030 2.54 
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0.083 2.50 
0.042 2.42 
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0.077 2.40 
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0.108 2.30 
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TABLE C-5 (continued) 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
130 0:035* 4.30 
0.062 4.25 
0.051% 4.10 
0.034 4.08 
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0.058 3296 
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TABLE C-5 (continued) 
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Axial Location (z/D) Radial Location (Yr) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
2 0.136 2.14 
0.123 2.14 
0.140 1.92 
0.159% 1.92 
0.138* 1.88 
0 134* 1/7 
0.147 1 6 
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Radial Location (r) 
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Radial Location (r) 
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Axial Velocity (Vz) 


(cm/sec) 
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TABLE C-6 (continued) 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
1.0 O.A137 5.41 
0.139 5.40 
O5125* 5.70 
0.037 5.20 
0.145* 4.95 
0.140* 4.80 
0.147% 4.35 
0.154* 6150 
0..165* 2.70 
ON i2% 2130 
2,0 0.000 7.95 
0.010 7.80 
0.000 7.80 
0.038 7.66 
0.013 7.60 
0.012 7.60 
0.040* 7.400 
0.045 7.42 
0.054 7.40 
0.029 7.40 
0.051 7.30 
0.056* 7.10 


0a. ¥: 
io. f 
fe. % 
sh. 
06,5 
oF - 
tie atk 


A Re | 
\  * ie 
veer Paes we ts. . 
irr “3 Te oe 


Axial Location (z/D) 
(No. of Dia.) 


2.0 


TABLE C-6 (continued) 


Radial Location (r) 
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Axial Velocity (Vz) 


(cm/sec) 
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TABLE C-6 (continued) 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
de O477* Cane 
0.189* 2705 
0% 196% 1.91 
3.0 0.000 5.86 
0.013 5.86 
0.000 i de he) 
0.042 5.62 
0.012 5.60 
0.048* 5.49 
0.061 5.47 
0.070* 5.34 
0.054 5720 
0.064 Shy 2 
0.015 ey 40) 
0.032 5.00 
07069" 5.00 
0.107 4.99 
0.074 4.95 
0.070* 4.78 
0.073 4.70 
0.121 4.66 
0.094 4.64 
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TABLE C-6 (continued) 


Axial Location (z/D) Radial Location (r) Axial Velocity (Vz) 
(No. of Dia.) (cm) (cm/sec) 
Sa) Oetty* 4.53 
0.100* 4.46 
OB 4.28 
Galose 4.20 
Oc i225 4.00 
0.144% Sie e 
0.138 3.80 
Ueto Broce 
0.170 3730 
0.1/8" 2.9] 
05713" 2.6/ 
Oise 2.66 
0.194% 2.66 
0.190 2.60 
02:210* 2205 
022267 1.90 


O2232> 1.60 


i 


nee tee ee ee ee 5 a cg SN. — 


143 


Cz Centerline Velocities 

Contained in this section are the axial velocities at the 
flow centerline determined from the plots of the axial velocity profiles 
such as Figures (IV-1) and (IV-2) at several upstream axial locations 
and for various flow rates. For example, considering the 0.30 cm I.D. 
capillary tube at flow rate 0.6392 cm?/sec and at axial location of two 
diameters (z/D = 2), the data points are plotted in Figure (IV-2). By 
drawing an eye-fit solid line through the data points, the centerline 
velocity was estimated at y = 0 and was found to be 7.95 cm/sec. This 
is the value tabulated in Table (C-7). 

For the axial locations where a complete velocity profile 
was not measured, the centerline velocities were determined by the same 
procedure except the incomplete velocity profiles near the centerline 


were used instead of the complete profiles. 
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TABLE C-7 
Diameter Flow Rate Axial Location Centerline Velocity 
(z/D) (Vz, ¢) 
(cm) (cm?/sec) (Nos-of- Dig.) (cm/sec) 
0.20 0.1598 — 0.0 8.25 
0.5 7332 
T.0 6.38 
lS Sey 
2.0 4.85 
2.5 4.00 
6-0 3.60 
0.3196 0.0 . 6.30 
0.5 15.00 
1.0 13.80 
15 12.60 
2.0 T1253 
ese 10.40 
3.0 OM55 
0.6392 0.0 32230 
30.30 

v0 28.01 
2.0 24.40 
3.0 Z\PL0 
4.0 18.00 
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TABLE C-7 (continued) 


Diameter Flow Rate Axial Location Centerline Velocity 

(z/D) (Vz,¢) 

(cm) (cm?/sec) (No. of Dia.) (cm/sec) 

0.20 0.6392 5.0 15350 
6.0 13.40 
fhe’ 11.55 

0.30 0.1598 0.0 3.60 
0.5 apt ayes 
beO 1585 
es yas 
2.0 0.85 

0.3196 0.0 fbeas 
005 6.0] 
1.0 4.75 
125 Ths. 
2.0 3.00 
2o5 2.26 
3.0 1.80 
0.6392 0.0 14.50 

0.5 12555 
10 10.82 
15 9.15 
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TABLE C-7 (continued) 


Flow Rate Axial Location Centerline Velocity 
(z/D) (Vz.¢) 
(cm? /sec) (No. of Dia.) (cm/sec) 
0.6392 220 795 
205 O57 5 


3.0 5.86 
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APPENDIX D 


CONE ANGLES AND CENTRAL CORE DIAMETERS 


The cone semi-angle was defined in Chapter III, Section 
III.3-3b, as the angle formed by the line tangent to the outermost 
streamline that enters into the capillary tube. The cone semi-angles 
which weremeasured in degrees and the central core diameters which were 
measured in centermeters are tabulated in Table (D-1) for the 0.20 and 
0.30 cm 1.D. capillary tubes at several axial locations upstream of the 
contraction. The initial cone semi-angle (9) was defined as the cone 


semi-angle at the contraction (z/D = 0). 


TABLE D-1 
Diameter Flow Rate Axial Location Central Core Cone 
(z/D) Diamter (d) Semi-Angle(¢) 
(cm) (cm?/sec) (No. of Dia.) (cm) (degrees) 
0.20 0.1598 0.0 0.200 O70 
Tz0 0.220 4.10 
2.0 0.240 4.40 
osu 0.276 Saou 
4.0 0.316 6.60 
el 6 0.364 9.70 
6.0 0.434 12.60 
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TABLE D-1 (continued) 


Diamter Flow Rate Axial Location Central Core Cone 
(z/D) Diameter (d) Semi-Angle (¢) 
(cm) (cm?/sec) (No. of Dia.) (cm) (degrees ) 
0.20 0.3196 0.0 0.200 2.30 
130 Os715 (dy toll) 
220 0.228 2.60 
320. 0.239 2.80 
4.0 0.261 Se 20 
0) 0.28] Sof AO) 
6.0 02305 4.30 
7.0 Oe3o/ bao 
0.6392 0.0 0.200 60 
1.0 0.209 1 66 
Zag Oris ‘bes! 
i.) 0.230 1.84 
4.0 0.24] 2210 
5-0 F205 2.40 
6.0 0.268 2.80 
720 0.287 Syed) 
0.30 0.1598 0.0 0.300 9.40 
O25 363 11.80 
1.0 0.432 15340 
ets) 0.526 19.40 
2.0 0.657 26.00 
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TABLE D-1 (continued) 


Diameter Flow Rate Axial Location Central Core Cone 

(z/D) Diameter (d) Semi-Angle (4) 
(cm) (cm?/sec) (No. of Dia.) (cm) (degrees) 

0.30 0.3196 0.0 0.300: 6.20 

O55 O2303 4.90 

FG ORS 72 8.00 

eo 0.416 9.20 

eat) 0.472 11.80 

20 0.546 ford 

a0 0.642 22.90 

0.6392 0.0 0.300 4.10 

O25 07325 4.50 

120 0.344 4.70 

15 0.369 4.90 

2.0 Ges9c 5.60 

250) 0.424 6.60 

SAU0) 0.464 7200 

Se) O2503 9.40 

4.0 0.560 ainso 


APPENDIX E 
STRETCH RATES ALONG CENTERLINE 


The stretch rates along the centerline were measured as the 
slope of the plot of axial centerline velocities (Vz,¢) versus the 
axial locations (z). The data which were measured are tabulated in 


Table (E-1) at several axial locations upstream of the contraction. 
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TABLE E-1 

Diameter Flow Rate Axial Location Centerline 
(z/D) Stretch Rate 

(cm) (cm?/sec) (No. of Dia.) (sec!) 

OF20 0.1598 035 = 10.39 

fe - 9.16 

1.5 - 8.12 

“aii = Pals 

2.9 - 6.06 

3.0 - 5.14 

0.3196 0.5 => o.00 

1.0 are | ro" 

105 - 11.45 

2.0 - 10.81 

feo - 9.78 

oa0 =1 0,99 


a a) 


Poneto 


ont; ghe7ay oan salted a 
tahatgtude!, ear tein ¢ 
Ast tasys nes gtd Ye wiartegy 


iN 


ee ee 


ra } 3 a) 


oo 


Fale 
a 


ea 


p>! | 
Cr. 


G.{ 


cH, 
— 


G4 
oy2 


oe 


By | 
TABLE E-1 (continued) 


Diameter Flow Rate Axial Location Centerline 
(z/D) Stretch Rate 
(cm) (cm?/sec) (No. of Dia.) (sec!) 
0.20 0.6392 075 - 21.71 
18 - 20.63 
(ah - 17.96 
3.0 - 15.64 
4.0 - 13.61 
one - 11.60 
6.0 - 10.00 
20 = Upollé 
0.30 0.1598 Uap = Pash) 
1.0 - 4,43 
ee a3 %\ 1 
2.10 = 2.26 
0.3196 25 =a el 
1. eT ys 
es - 6.05 
Ae) - 4,94 
Ceo macs 07, 
CAL! - 2.60 
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TABLE E-1 (continued) 


Diameter Flow Rate Axial Location Centerline 
(z/D) Stretch Rate 
(cm) (cm?/sec) (No. of Dia.) (sec!) 
0.30 0.6392 0.5 - 12.84 
1.0 - 11.24 
1.5 - 9.6] 
2-0 - 8.19 
25 - 6.86 
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APPENDIX F 


SAMPLE CALCULATION 


Contained in this section is a sampie calculation which converts 
the measured streak length into axial velocity. For example, consider 


the following: 


TD. *ofcapiitary Stube’vs=4 00.30% om 
flow rate = 0.3196 cm? /sec 
rotational speed of the 20 slits and spokes chopping disc 
=a R Pee 
magnification of streaks project onto the graph paper 


=e 


Typical streak lengths along several streamlines for this run 
were shown in Figure (III-6). Consider the streamline number 8 in 
Figure (III-6) at the axial location of one diameter (z/D = 1). 

(i) The streak lenths were measured by using the leading edge 


to leading edge technique or distance A to B. 


In this case the length of AB was equal to 3.70 cm on the 
graph paper. Then the length of streak (AL) at this axial 


location can be found by dividing this value by 2. 


AL = TerTRS i 1.85 (cm) 
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(ii) The time interval for this example can be calculated by 


using the following Equation (G-1) 


fe at OUR eEe (G-1) 
Nw 
where 
N = number of slit and spokes on the chopping disc 
w = number of revolution per minute. 


Therefore, 


60 


t 70 x 242 


0.0124 sec 


(iii) By using Equation (III-1), the axial velocity can be 


calculated. 

Vaan (III-1) 

Mt 

where 

AL = length of streak 

M = magnification of streaks on the graph paper 

t = time for the slit to cut the beam of light 
Therefore, 

Vz = 185 


47.8 x 0.0124 


3.lc; ~em/sec 


This axial velocity was plotted in Figure (III-7) and also 


tabulated in Appendix C, Table (C-5). 
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